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We construct time-periodic solutions for a system of self-gravitating massless scalar field, with 
negative cosmological constant, in d + 1 spacetime dimensions at spherical symmetry, both pertur- 
batively and numerically. We estimate the convergence radius of the formally obtained perturbative 
series and argue that it is greater then zero. Moreover, this estimate coincides with the boundary 
of the convergence domain of our numerical method and the threshold for the black-hole formation. 
Then we confirm our results with a direct numerical evolution. This also gives strong evidence for 
nonlinear stability of the constructed time-periodic solutions. 



Introduction. A recent numerical and analytical study 
of the four dimensional spherically symmetric Einstein- 
massless scalar field equations with negative cosmolog- 
ical constant indicated that anti-de Sitter (AdS) space 
is unstable against the formation of a black hole under 
arbitrarily small generic perturbations pQ. Qualitatively 
the same results were obtained later in higher dimen- 
sions |2 3 . Although gravitational collapse seems to be 
a generic fate of a small perturbation of AdS, it was sug- 
gested in [T] that there may exist non-generic initial data 
for which the evolution remains globally regular in time. 
This conjecture was later generalized to the vacuum Ein- 
stein equations with negative cosmological constant [4j [5] , 
where time-periodic solutions (geons) were postulated. 
In this note, continuing the study initiated in pQ, we 
give strong evidence for the existence of time-periodic 
solutions of the spherically symmetric Einstein-massless 
scalar field equations with negative cosmological constant 
in d + 1 dimensions for any d > 2. We construct these 
solutions using two independent methods: nonlinear per- 
turbative expansion and fully nonlinear numerical evolu- 
tion. The fact that these two methods produce the same 
solutions make us feel confident in our results. In the case 
of even d the perturbative construction of time-periodic 
solution can be algorithmized with the highest order of 
the expansion limited only by the memory resources of 
a computer on which our Mathematica script is running. 
The case of odd d is more involved due to incompatibil- 
ity between the boundary behavior of the eigenmodes of 
the linearized problem and the solutions of the full sys- 
tem and its perturbative treatment will not be discussed 
here (time-periodic solution can be still constructed with 
slightly different method 0). On the other hand our nu- 
merical method for constructing time-periodic solutions 
works for any d > 2. 

After constructing the time-periodic solutions we check 
our results through direct numerical evolution of their 
initial data. The fact that this numerical evolution re- 
produces those time-periodic solutions serves as another 
cross-check of our results but, even more importantly, it 
is a strong evidence that these time-periodic solutions are 
(nonlinearly) stable. Thus, an interesting picture of dy- 
namics in asymptotically AdS space emerges: although 



AdS space is generically unstable against black hole for- 
mation, it possess also islands of stability. This picture 
was previously advocated in [5]. 

It is well known that in asymptotically flat spacetimes 
there are no non-trivial time-periodic solutions [7J |8]. 
The mechanism that allows for non-trivial time-periodic 
solution in our case is the lack of dissipation of energy. 
The existence of stable time-periodic solutions of the Ein- 
stein equations with negative cosmological constant is 
interesting in its own right, but it would be also fasci- 
nating to learn what is the counterpart of their stability 
islands on the gauge theory side, using AdS/CFT corre- 
spondence. 

The finite difference code used in [T] was well suited to 
perform stable, long-time evolution for the initial data 
that ultimately collapsed to a black hole, but in the case 
of solutions that stay smooth forever (like time-periodic 
solutions) spectral methods are more efficient and more- 
over they also allow for a direct comparison with pertur- 
bative construction of time-periodic solutions. Thus, we 
construct a pseudo-spectral code that is well suited to 
deal with the problem at hand. 

Finally let us mention that the time-periodic solutions 
we construct form a particular class of time-periodic so- 
lutions with one dominant mode. Whether other time- 
periodic solutions, that do not bifurcate from one-mode 
solutions of the linearized problem, exist or not is an open 
question. 

In this note we deal with the simplest model case of 
the self-gravitating massless scalar field, but the meth- 
ods presented here can be also applied (with minor mod- 
ifications) to the pure vacuum case, at least for the 
cohomogeneity-two Bianchi IX ansatz in even number 
of spatial dimensions, discussed in [9]. 

Model. To make this note self-contained we rewrite 
the equations for the selgravitating massless scalar field 
with negative cosmological constant from [TJ We 
parametrize the (d+ l)-dimensional asymptotically AdS 
metric by the ansatz 

f 2 

ds 2 = 7T- (~Ae~ 2S dt 2 + A~ 1 dx 2 + sin 2 a; dtt 2 , ,) , (1) 

cos 2 a; v « ±/ ' w 

where £ 2 — —d(d — 1)/(2A), d£l 2 d l is the round metric 



2 



on S^ 1 , — oo < i < oo, < x < n/2, and A, S are 
functions of (t, x). For this ansatz the evolution of a self- 
gravitating massless scalar field 0(t, x) is governed by the 
following system (using units where 8irG = d — 1) 



$ = (Ae- s Il) 



n 



i 



tan a; 



The set of equations ( [l0| [5| [6j [9| forms a closed system 
in the form well suited for the numerical integration. 

Linear perturbations. The spectrum of the linear self- 
adjoint operator, which governs linearized perturbations 
, of AdSd+i, L = — (t&nx) 1 ~ d d x [(tanx) d ~ 1 d x ) , is given 

— (tan d_1 a; Ae~ $) , (2) by w 2 = (d + 2j) 2 , j = 0,1,.... The eigenfunctions read 



A'-- 



d-2 + 2sin z x 



sin x cos x 



5' = - sin x cos x ($ 2 + II 2 ) , 
where ' = dt, ' = d x , and 
3> = </>', 



(1 - A) - sin a; cos x A ($ 2 + IT 2 ) , 

(3) 
(4) 



n = A"Vc 



(5) 



The key point is to note that in this particular gauge both 
constraint equations ([3j [4| can be put in the form that 
allows for efficient integration. This idea was previously 
exploited in |10) for the efficient parallelization in solving 
the constraints in finite difference code, but here we will 
use this form to construct a spectral code and to find the 
time-periodic solutions of the system both perurbatively 
and numerically. For a given matter content, prescribed 
with the functions and IT treated as independent dy- 
namical variables, the metric function S can be expressed 
as the integral 



3in y cos y ($ 2 [t, y) + IT 2 {t, y)) dy , (6) 



where the value of xo reflects the residual gauge choice 
in ([!]) . The choice xq = makes the time coordinate t to 
be the proper time of the central observer. Now, using 
Q and the identity 



d - 2 + 2 sin x cos d x ( sin x 



sin x cos x 



sin d 2 x\ cos d x 



we can rewrite ([3| as 



cos" x \ cos d x j \ cos d x 



(7) 



(8) 



Multiplying this equation by e s and integrating by parts 
we get 



1 — A = e 



(cosx) d 



(sinx) d 2 

x / e- s ^^ 2 (t, y )+n 2 (t,y))(t^y) d - 1 dy. 



(9) 



It is also convenient to rewrite the wave equation ([2| 
eliminating A' and 6' in the equation for IT, using the 
constraint equations (pi El) . In this way we get 



w t / , ^ d - 1 - (1 - A) cos 2x\ 
IT = e d \&A + $ ^ ). (10) 



sin x cos x 



,(x) = 2 ^!0- + ^l)! (cosx) ^/2-M/2 ( 



cos 2x) , 



(11) 

where P"'^(x) are the Jacobi polynomials. These eigen- 
functions form an orthonormal base in the Hilbert space 
of functions L 2 ([0,7r/2], (tanx) 4 " 1 dx) . Below we de- 
note the inner product on this Hilbert space by (/ \g) := 
J^ 2 f(x)g(x)(tanx) d ^ 1 dx. It will be useful in the fol- 
lowing to note that (e'^x) \e'j(x)) = u 2 8ij. 

Perturbative construction of time-periodic solution. We 
seek a time-periodic solution of the system (|2]|4| in the 
form 

cj>= e cos(r)e 7 (a;) + e x 4>\(t, x), (12) 

odd A>3 

5= J2 eA ^(r,x), l-A= J2 £X Mr,x), 



even A>2 



even A>2 



(13) 



where e 7 (x) is a dominant mode in the solution in the 
limit e —¥ 0, r = Cl^t is the rescaled time variable with 



f2 7 = cj 7 + /J 



W 7 X 



(14) 



even A>2 



and 



4>x=Y,fxAr)ej{x), (15) 
j 

S\ = <Za,-i(t) + 22dx,j(r)ej(x), A x = }^ax,j(T)ej(x) . 
j j 

(16) 

with f\,j(r), a\j{r), d\j(r) being periodic in r. It is im- 



portant to note that for even d the sums in (15 16 ) are 



finite at each order A of the perturbative expansions ( 12 
13 1. This allows for a straightforward algorithmization of 



building up the successive terms in (12 13 1. That is not 



the case for odd d. This is due to the boundary expan- 
sion of the solutions of the system (2|4|. Smoothness at 
spatial infinity and finiteness of the total mass M imply 
that near x — n/2 we must have (using p = n/2 — x) 



4>{t,x)=0{p d ), 

A(t,x)=l + (D(p d ), 6(t,x) 



0( P 



2d\ 



(17) 



It can be checked that for odd d, higher even terms in 
the expansion of <fi at the boundary x = n/2 do not 
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vanish [TT] and cj>\ does not have finite expansion into 
eigenmodes ej(x), being (in odd d) odd, with respect to 
x = tt/2. For later convenience it is useful to introduce 
the following notation: let [e A ] / denote the coefficient at 
e x in the (formal) power series expansion of / = J2x 



Now, inserting the first of the series in (13 16 ) into Q 
and projecting onto e' k {x) we get 



dx.k = - 



2^ 



[e x ] sin(2a;) ($ 2 + n 2 ) ) 



(18) 



and cI\—i(t) is fixed by the gauge fixing condition 
6\(t, 0) = 0. Similarly, inserting the last of the series 
m ([13] [16]) into ([8f using Q [7] and then projecting onto 
efc(x) we get a linear system of equations for the coeffi- 
cients a\ t j{r): 



E 



(d - l)S kj + e k 



— sin 2x e'j — cos 2x e.j 



"A, 



(e k \[e x ] (sin2a;) 2 A($ 2 + n 2 



(19) 



It is useful to note that the principal matrix of this 
system is 3-diagonal. This system supplied with 
[e A ] (1 — A)\ _ Q = = J2j a Aj' e i(0) condition allows for 
the unique solution for the coefficients a\j(r). Then, 
for odd A > 3, one gets that (f>\ fulfills an inhomo- 
geneous wave equation on the pure AdS background: 
(w 2 <9 TT + L) 4>\ = S\. Projecting this equation onto 



efc one finds that the coefficients fx.k in (|15 ) behave 



as forced harmonic oscillators: [u>^d TT + us k 
(ek\S\). Solving these, we get two integration con- 
stants for each of the equations. The two constants 
in /x ~ (the coefficient for the dominant mode e 7 in 



(15l), are fixed with the condition (/a. 7 , d T fx,~f)\ T=0 



(070), that is we choose ((e 7 \<f>) , (e 7 \9 T (p))\ T _ 
(e, 0). Now, if (ek\S\) contains the resonant terms 
cos(wfc/w 7 )r or sin(wfc/w 7 )r, this gives rise to secular 
terms rsin(rwfc/w 7 ) and TCOs(rajfc/w 7 ) in fx .k respec- 
tively. Such terms would spoil the periodicity of the 
solution, thus they have to be removed. This fixes 
the correction to the frequency w 7! a-i, and the inte- 
gration constants. Namely, it turns out that in order 
to not to produce spurious resonant terms in higher or- 
ders, all but odd (in r) frequencies in the solutions for 
fx,k have to be removed and moreover, all fx,k tune 
in phase to the dominant mode: d T fx,k\ T=0 = (for 
the choice <9 r /A l7 | r=0 = 0). In summary, at any odd 
A > 3 we fix w 7 .a-i and, as (e^ \Sx) does not vanish 
only for a finite number of modes (namely (e^ \Sx) = 
for k > 7 + (d + 1 + 27) (A - l)/2), we are left with 
(A - l)/2 + L( A - 1 )/( 2 (rf + 2 7))J undetermined integra- 
tion constants. They will be fixed, together with w 7 .a+i, 
to remove (A + l)/2+ [(A- l)/(2(d + 27))J secular terms 
in 0a+2- 

In fact all the projections onto (or e' k ) appear- 
ing at any order of the perturbative procedure de- 



scribed above, can be reduced to just a few in- 
ner products: (efe|e,ej), (ek |cos 2x Ci ), (e^ |sin 2x e\ ), 
(ejt |csc x sec x cos 2x e'j ) , (e/. |csca; secxei e'j). Thus, 
the whole procedure of building up such perturbative so- 
lution is straightforward to implement. 

Pseudo-spectral code for the time evolution. We expand 



and II into K eigenmodes of linearized problem (11) 



E > 

0<j<K 



fj(t)ej(x). 



n 



P 3 (t)e 3 (x) (20) 

0<j<K 



Prescribing the coefficients fj(t), Pj(t) at time t, we want 
to solve for them at a later time t + dt. To achieve that, 
we calculate their time derivatives fj(t) and pj(t) with 
a spectral method and then integrate them in time. We 
know that the metric function S and the integrand in ^ 
can be efficiently [12] expanded as follows 



-8 



5 = \_. dj(t) cos(2jx) . 

0<j<K 

($ 2 + n 2 )= ^2 Uj (t) cos(2jx) . 



(21) 
(22) 

n 2 ) . 

(23) 

When the sides of ( 23 1 are evaluated at K collocation 
points, Xk — (ir/2)k/(K + 1), k = 1, . . . , K, we get the 
linear system of equations to be solved for the coefficients 
dj(t). As da(t) is absent in this system, reflecting the 
gauge freedom, we add the gauge fixing condition 



Q<j<K 

Substituting ([2l]) into Q we get 

(— 2j) sm(2jx)dj(t) = — sinx cosx ($ 2 

l<j<K 



E "/ 

0<j<K 



d 3 (t) = 0. 



(24) 



Then, evaluating the sides of (22 1 at the collocation 



points we get the linear system of equations to be solved 
for the coefficients Uj(t). This allows us to solve for the 
metric function A. Using |9]) we get 



A = 1 



Si E «*W*> 



(sinx)'*- 2 ^ 
v ' 0<j<K 



(25) 



where the weight functions (x) read: 



w i ( x ) = J cos ( 2 i2/)( tan y) 



d-l 



dy. 



(26) 



Now, substituting the expansions ( 20 ) into wave equation 
( 10 5]) and evaluating both sides at the collocation points 
we get the linear system of equations to be solved for the 
time derrivatives fj(t) and pj(t). 

Numerical construction of time-periodic solutions. Seek- 
ing for time-periodic solutions numerically it is conve- 
nient to use rescaled time coordinate t = fit where, as in 



4 



the perturbative construction, f2 denotes the frequency 
of the solution we are looking for. Assuming that time- 
periodic solution does exist, we expand <p(t, x) and IT(i, x) 
into eigenmodes of the linearized problem in space and 
Fourier coefficients in time. Choosing a grid with K col- 
location points Xk = (n/2)k/(K + 1), k = 1, . . . , K in 
space (as for the spectral code) and N collocation points 



in time r„ = 7r(n - 1/2)/ (2JV + 1), n 
truncate these expansions as follows 



1.. 



<t>= kj cos (( 2i + l ) T ) e o( x ) 

Q<i<N 0<j<K 

n= Yl E ft,iSin((2i+l)r)e J -( a: ) 

0<i<JV 0<j<K 



,N we 

(27) 
(28) 



Next, at each instant of time r n we calculate the coeffi- 
cients 



fj(tn)= Yl cos((2i + 1)t„) , (29) 

0<i<JV 

&•(*») = 1] p iJ sin((2i + l)r„) (30) 



0<i<N 



and put them as an input to our spectral procedure, 
getting as the output their time derivatives fj(t n ) and 
Pj(t n ). Equating those to the time derrivatives of (28) 



(remembering that dt = £1<9 T ) at the set of K x N grid 
points (r„, Xk), n — 1, . . . , N, k = 1, . . . ,K, together 
with (28) inserted into the second of the equations ^ 
and then evaluated at the set of K x N grid points, 
and the equation ^2 0<i<N /i, 7 = e, setting the ampli- 
tude of the dominant mode 7 in the initial data to e, we 
get a nonlinear system of 2x A"x JV + 1 equations for 
2 x K x N + 1 unknowns: f+ j, Pij and 57. This system 
is solved with Newton-Raphson algorithm yielding the 
time-periodic solution for the system (2]|4). 
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8.72677- 10" 
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-1.65282 ■ 


10" 
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9 


2.74313 ■ 10" 


14 


63 


2.31645 ■ 10" 


69 







TABLE I: Coefficients of the expansion ( 20 \ at t 
(with Pj(0) = 0), corresponding to the time-periodic solu- 
tion with dominant eigenmode 7 = 0, with the amplitude 
e = 1/100, determined by the numerical procedure with 
K x N = 64 x 24 grid points and extended floating-point 
precision. The solution oscillates with frequency ~ 4.00667. 



Results. The detailed results obtained with the tools de- 
scribed above will be given in a longer accompanying 
paper. Here, just to give an example, we present the 



time-periodic solution bifurcating from the fundamental 
mode (as a dominant one) in 4 + 1 dimensions. First, 
we construct the time-periodic solution with a numerical 







FIG. 1: Visualization of the phase space generated by time 
evolution of initial data listed in Tab. [I] All plots depict the 
same solution evolved over time interval equal to 500 peri- 
ods. Closed curves on the slices of phase space give a strong 
argument for the stability of analyzed solution. 



procedure described above. Then, we read off the coeffi- 
cients of the expansion (20) at t — (see Tab.[l]) and put 
them as the initial data into our spectral evolution code. 
Despite the presence of truncation errors and numerical 
noise in these initial data, their time evolution is peri- 
odic in time, as depicted by closed loops in Fig. [T] repre- 
senting different sections of the phase space, spanned by 
the set of coefficients {fj(t), Pk(t)}- This provides strong 
evidence not only for the existence of the time-periodic 
solutions but also for their (nonlinear) stability. This ar- 
gument for the stability is strengthen by the fact that 
if we perturb the time-periodic solution slightly then its 
evolution is no longer periodic, but it does not collapse 
to a black hole and stays close to the periodic orbit. Us- 
ing the perturbative method, for small values of e we get 
fast convergence to the numerical solution, for larger e 
the convergence can be improved with the Pade resum- 
mation (see Tab. [Hj. The Pade approximation can be 
also used to estimate the convergence radius for pertur- 
bative expansions ( 12pT ). For example, if we construct 
the (n, n) Pade approximants for 
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^ 464 2 45614896 4 173158711507904383595696 6 19627018631453126466665156076805265104 8 
+ 6 + 11319 £ + 533797475350414275 £ + 662148921092395909349993941125 

11072083972904297030696081837311640731963665649719830983726203187072 10 
+ 3798775287765327851156088282697847314054859440405429239375 £ 

20613105125948158448048760017030229079689748556646392089437804474799052885011038708688 12 
+ 68546459749514836602259551946772532940114419607536660561179461719235846875 £ 

+ uj 0M e u + 0J ,i6S 16 + O(e 18 ), 

(31) 



with n = 2, 4, 6, 8 [13] and then calculate the zeros of the 
denominators nearest to the origin we get 0.128, 0.102, 
0.095 and 0.092, respectively. These values are consistent 
with the fact that the Newton-Raphson algorithm used 
in the numerical construction of time-periodic solutions 
ceases to converge for e > 0.09. Moreover if we read off 
the initial data for larger values of e the time evolution 
leads to a black hole formation almost immediately. 



Padc 



0.005 
0.015 
0.025 
0.035 
0.045 
0.055 
0.065 
0.075 
0.085 



4.0016596666501 

4.0151220741462 

4.0430867838460 

4.0879197007 

4.15407139 

4.249920 

4.39267 

4.6230 

5.05 



4.0016596666501 

4.0151220741462 

4.0430867838521 

4.0879197035435 

4.15407167953 

4.249932516 

4.3929928 

4.629225 

5.18 



4.0016596666501 
4.0151220741462 
4.0430867838521 
4.0879197035448 
4.1540716797440 
4.2499325336279 
4.3929938556099 
4.6292962269712 
5.2017714694183 



TABLE II: Frequency of solutions with 7 = determined 
by: direct summation of the perturbative expansion fis | |31| , 
diagonal (8,8) Pade resumation f2p a dc and the numerical code 
!] nllm . The relative absolute difference of Q nU m and fip a de 
ranges from 10" 28 for e = 0.005 to 4 ■ 10" 3 for e = 0.085. 
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